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Leakage current is one of the critical aspects to consider for overhead 
transmission line insulator’s condition and performance assessment. 
As the leakage current increase, the size of the dry band also will increase leads 
to the dry band arcing, deteriorate the insulator performance and contribute to 
the development of insulator flashover. Based on the literature study, other than 
the existence of contaminations on the surface of insulator combined with 
moisture, the variation in leakage current is also affected by the environmental 
and electrical stresses. Previous researches have shown the effect of 
environmental and electrical stresses on surface leakage current based on 
experimental and simulation results. This paper outlines an analytical approach 
based on dimensional analysis to propose a new mathematical model of 





Leakage current (LC) leakage current under environmental and electrical stresses. To justify 
Mathematical modeling the applicability of the derived dimensional model, the new model has been 
contamination validated using previous researcher’s experimental results. The validation 
indicated that the proposed model had shown a good agreement with 
the previous experimental results. The proposed dimensional equation for this 
research work can be potentially used as a predictive performance model to 
evaluate and monitor the leakage current and insulator’s performance. 
This is an open access article under the CC BY-SA license. 
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1. INTRODUCTION 

Transmission line insulator, also known as overhead insulator is greatly affected by the environmental 
stresses produced by the weather and surroundings. Insulator is usually exposed to uncertain weather 
and environmental pollution [1] such as inert material, soluble salt, water, electronic conductive dust, 
chemicals, industrial contaminants, and agricultural pollutants [2, 3]. Due to this variety of pollutions, 
insulator’s external surface becomes contaminated and tends to accumulate the contamination, allowing 
the leakage current (LC) to flow. The contamination deposited on the insulator surface combined with 
the presence of moisture dissolves the salt deposited and encourages the formation of the conductive layer 
that acts as a highly variable and nonlinear resistor to be formed for LC to flow [3, 4]. With high LC flowing 
between the high voltage and grounded side of insulator, flashover occurrence that will lead to power outage 
from the system is more likely to happen [5-7]. Flashover phenomenon reduces the power system reliability 
and efficiency, contributes to economic losses, and permanent insulator impairment [8, 9]. Due to these 
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reasons, LC monitoring is one of convincing way to examine the insulator performance especially under 
contaminated conditions [2, 10, 11]. LC monitoring describes the insulator’s surface condition and is capable 
to act as an accurate indication tool for flashover warning and hence prolonging the insulator’s lifetime [12, 13]. 

In recent years, researchers have investigated the effect of environmental and electrical factors on 
LC. However, most of the research results come from experiments and simulations. Testing on high voltage 
prototype sometimes causes costly mistakes, requires a relatively long considerable time, and there are cases 
when the analytic expression of the variables is not available or inaccurately known [14, 15]. Some research 
efforts were also directed towards establishing a mathematical correlation and modeling, which is adapted 
from empirical approach [16-20]. However, the models developed only effective as the LC prediction 
and unable to represent the actual physical equation that explains the characteristics of the LC. This paper 
presents the mathematical relationship to describe the LC, considering the environmental and electrical 
factors using the analytical approach of dimensional analysis (DA). DA is an analytical technique which 
applies the concept of similarity to explain the relationship among the physical quantities [14, 15, 21]. 
By using DA approach, the physical phenomenon is represented by dimensionally correct equation among 
certain variables [22, 23]. 


2. DEVELOPMENT OF DIMENSIONAL MODEL 

From the literature studies, it can be concluded that the LC characteristics have a strong relation 
with many physical variables from the environmental and electrical stresses. The LC that flows on insulator’s 
surface depends on the dominant variables such as Equivalent Salt Deposits Density (ESDD), humidity, 
electrolyte conductivity, environmental temperature, applied voltage, environmental pressure, ultraviolet 
radiation, and insulator’s creepage distance. A mathematical relationship to describe the LC is established 
between these variables using DA. The base dimension of the target variables is shown in Table |. Among 
the target variables, I, is considered as a dependent variable and Esp, h, Og, Ta, Vs, Pe, Suy and Lç as 
independent variables. The unit system used in developing the LC mathematical model is SI unit standard 
system and the base dimensions are mass (M), length (L), time (T), charge (Q) and temperature (0). 


Table 1. Base dimension of the target variables 








Variables Symbol Base dimension 
Leakage current (LC) I, QT 
Equivalent salt deposits density (ESDD) Esp ML? 

Humidity h ML 
Electrolyte conductivity Og M*L3TQ? 
Environmental temperature Ta 0 

Applied voltage Vs MLT Q 
Environmental pressure Pe MLT? 
Ultraviolet radiation Sup MT 

Insulator creepage distance Tig M*L"T2Q 





From the target variables, the relationship between the independent and dependent variables are established 
as follows, where f is an unknown function. 


I, = f(Esp.h, Or, Ta, Vs, Pes Suv» Lc) (d) 


From the list of the target variables, the number of variables, Ny = 9 and the number of dimensions, 
Na = 5. DA adopts the Pi Buckingham’s theorem in its analysis. Buckingham method starts by examining 
the number of variables involved and the number of dimensions [15]. The final result is a relation between 
Ny and Na which provides the number of dimensionless variables (no of II’s). Therefore, by applying Pi 
Buckingham’s Theorem, the no of independent set of product (IT’s) is IP’s= Ny — Na = 9 — 5 = 4. The 
dimensionless product can be written in a form as in (2). 


IL = Eçp®t. h2. opê3. Tat. Vs. Pe. Suy. Le“ 2) 


where €j aea €g are the exponent of the corresponding target variables. The dimensional matrix of 
the variables correspond to their base dimensions with MLTQ@ combination is arranged as follow; 
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The rank of the dimensional matrix is established as 5. To compute the independent set of product, 
the dimensional matrix is partitioned into two submatrices, Matrix A and Matrix B. Matrix A must be a square 
matrix and not singular. If Matrix A is singular, there is no possibility to produce the independent set of 
products (II’s). The number of independent set of product (II’s) is generated by using the formula (3)-(4) [15]. 


P=E.Z (3) 
z=] and He eis ges (4) 


where: P_ : Independent set of product 
B : Matrix B 

A : Matrix A 

q : Exponent of dimension 

I : Identity matrix 

e€ : Exponent of variable 

0 : Null matrix 


Matrix P established in (5) described the target variable and product of variables [15]. 


1 0 0o 0 
0 1 0 0 
0 oO 1 0 
0 oO 0 1 
Se ae ee a (5) 
Sa A PY i 1 
4% -% % i 
Df IBS 28. 29, 
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Tee Te ON 
The last step of the DA modeling is to produce the set of dimensionless target variables. The 
complete set of dimensionless products is now rewritten as in (6)-(9). 


2 





_ IL-Pe5 
Tı =- 7 (6) 
Vẹ-Suv-L? 
t 
= (Esp) Suv? Le 7 
2S. <4 aa (7) 
Ve.Pe 5 
hSyy? 
mt, = (Su (8) 
Pe 
r 
5 
DEV 
T4 =— a (9) 
Pe5-Suv-LÈ 


By referring to the Buckingham’s Theorem, the established dimensionless product is given as in (10). 
I = f Mh H, H, I4) (10) 


To show the relationship between the LC and other independent variables, monomial power 
form [24] is applied as shown in (11). Factors œt, œ? and a? is assigned by considering the relationship 
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between the dominant variables and leakage current. It is found that the leakage current is directly 
proportional to ESDD, relative humidity, and applied voltage. 


1, = f (12%. 13%. 1,%?) where a, = 1,a, = 1 anda, =1 (11) 


Equation (11) is rearranged to show the relationship of LC to the other dominant variables. 


1 
Vs Suy?-Lo*\5 
p = (E). f. n.m) (12) 

Finally, the resulting model of LC characteristics can be expressed as a function of environmental 
physical variables in (13) where Dç is the dimensional constant and can be defined from the experiment. 
The dimensional constant established in this research will vary depends on experimental conditions. 


1 
y 1 
I, = De(Esp-hħ. og- Suy“) ri (13) 

From (13), it can be proved that the derived dimensional equation obey the dimensional homogeneity 
principle. The principle of dimensional homogeneity helps in describing the consistency and completeness of 
an equation in physical algebra [21]. Therefore, (13) is projected to be used as a predictive model in 
estimating the LC in this research. 


3. VALIDATION OF THE PROPOSED DIMENSIONAL MODEL 

For validation purposes, the proposed leakage current dimensional model is verified with other 
researcher’s experimental results. The value of the dimensional constant or test parameter is recommended to 
fit the proposed model with other researcher’s experimental results by using Mathcad Software. This section 
discusses the verification of the model at different ESDD, electrolyte conductivity, relative humidity, applied 
voltage and insulator’s creepage distance. 


3.1. Validation of the model at different electrolyte conductivity 

Validation of the proposed model at different electrolyte conductivity is achieved by comparing it to 
the experimental results, Suwarno et al. [25] on polymeric insulator subjected to inclined plane test under 
various levels of kaolin-salt pollutions. All the dominants variables (Esp, ħ, Vs, De, Syy, Lc) in (13) are 
assumed to be constant except the levels of electrolyte conductivity. By assuming all the mentioned variables 
are constant, the proposed model can be simplified as follows; 


I,(o) = D,.0 +1, (14) 


where D, is the recommended dimensional constant to fit in the model with experimental results. Equation (14) 
is altered by adding J, to adapt the experimental result. Based on the simplified model, it is found that 
the LC flows on insulator surface is significantly increase with electrolyte conductivity. Figure 1 shows 
the comparison of the LC from the model developed and the doted data from the experimental result [25]. 
It can be observed that the proposed model fitting is very close to the experimental results. Table 2 shows 
the proposed LC equation in relation to the electrolyte conductivity with its data error and standard deviation. 



































5 r r r Table 2. Empirical equation of LC in relation with 
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Figure 1. Comparison of LC of experimental 
and proposed model at different electrolyte 
conductivity [25] 
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From Table 2, it shows that the model has a quite low deviation of 5.92% with reference to 
experimental surface LC. By varying D,, the curve can be adjusted to suit the experimental data. This 
provides an approximation to form the equation between LC and ESDD as in Table 2. 


3.2. Validation of the model at different ESDD and applied voltage 

The proposed model is verified with the experimental data, Banik et al. [26], at different ESDD 
and applied voltage. Experiment was performed on porcelain insulator to measure surface LC with different 
pollutions conductivity and increasing of applied voltage in steps of 5 kV to maximum 40 kV under control 
environmental [26]. The independent physical variables except for ESDD and applied voltage are assumed to 
be constant. Hence, (13) can be formulated in a function of ESDD and applied voltage. 


1 
L (Esp, Vs) = D2. Esp. Vs5 + Io (15) 


By referring Figure 2, it can be observed that the model is in complete agreement with the previous 
experimental results [26]. 
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Figure 2. Comparison between model and experimental result of LC at different ESDD 
and applied voltage [26], (a) 10 kV, (b) 20 kV 


The model able to compute the LC for different ESDD and applied voltage with deviation of 6.04% 
and 11.86% for 10 kV and 20 kV of applied voltage respectively (refer Table 3). The LC is in linear 
relationship with ESDD and applied voltage. Increasing the applied voltage will increase the electric field 
distribution and affect the LC variation [27, 28]. The dimensional constant of D, also shows an increment in 
its value with increasing of the applied voltage. 


Table 3. Empirical equation of LC in relation with ESDD and applied voltage 


iri ; £ 1 
Empirical equation l, = 0.57Esp. Vs5 + 0.18 1, = 3.15(Esp). V5 + 0.28 


(Applied voltage 10 kV) (Applied voltage 20 kV) 





Error 0.0604 0.1186 
Standard deviation (Model) 0.0956 0.6066 
Standard deviation (Exp) 0.1118 0.7127 





3.3. Validation of the model at different relative humidity and ESDD 

Equation (13) is written in a function of relative humidity and ESDD as shown below. From (16), it 
is found that LC is in linear proportion with relative humidity and ESDD. LC increases with the increasing of 
contamination severity and relative humidity. Increasing of contamination severity combine with high 
relative humidity add up the moisture content, reduce the insulator surface resistance and lead to a high LC 
flows on insulator surface [29]. 


IL, (h, Esp) = D3. Esp. h+I1, (16) 


The proposed model is verified with experimental results, Tousi et al. [2]. Figure 3 shows 
the verification of the proposed model using the data from the experiment [2]. It can be observed that 
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the deviation of the model is quite low (see Table 4) which are 4.51% and 4.01% compared to the experimental 
result. This identify that the model is capable to predict the LC at different relative humidity and ESDD. 
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Figure 3. Verification of proposed model at different relative humidity and ESDD using the results from 
experimental [2], (a) ESDD 0.26 mg/cm’, (b) ESDD 0.358 mg/cm? 


Table 4. Empirical equation of LC in relation with relative humidity and ESDD 








Empirical equation l, = 1.71Esp.h + 0.12 I, = 1.36Esp. h + 0.16 
ESDD 0.26 mg/cm? ESDD 0.358 mg/cm? 

Error 0.0451 0.0401 

Standard deviation (Model) 0.0889 0.0974 

Standard deviation (Exp) 0.1102 0.1124 





3.4. Validation of the model at different insulator’s creepage distance 
To validate the LC model at different insulator’s creepage distance, (13) is simplified in a function 


of specific creepage distance as shown in (17) by assuming all independent variables except creepage 
distance as constant. 


6 
I, (Lc) = D4- LÈ + l, (17) 


The model is compared to previous research data, Parihar et al. [30], which presents the leakage 
current on composite long rod insulator with different specific creepage distance and contamination severity. 
Figure 4 shows the LC in dependency on the specific creepage distance of the proposed model and 
the previous research data [30]. From the graph, it can be observed that the LC is inversely proportion to 
the specific creepage distance. With a higher value of specific creepage distance, conducting layer could be 
prevented thus decreases the LC flows on insulator surface. From Table 5, it may simply have verified that 
the model has a very close agreement with the previous experimental results [30] with a deviation of 0.23%. 
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Figure 4. Comparison of model and previous research 
result at different insulator’s creepage distance [30] 
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4. CONCLUSION 

The present article describes an analytical approach on modeling the LC of insulator under 
environmental and electrical stresses based on DA technique. The proposed dimensional model derived in 
this article complies with homogeneity condition and expressible by a dimensionally homogeneous equation 
in terms of specified parameters. The proposed method also able to present the dominant parameters 
associated to the LC behavior and recommends the appropriate correlation of the test parameter value to be 
used when executing the test. Throughout this paper, it has been observed that the dimensional model has 
a good agreement with the experimental results. This clearly shows that the model able to estimate 
the surface leakage current on insulators under the influence of environmental and electrical stresses 
and contribute to an indication mechanism of flashover warning. Furthermore, the proposed work also might 
help in reducing the cost of testing and labor work for the study of the insulator subjected to environmental 
and electrical factors. 
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